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It is shown that minimally coupled scalar field in Brans-Dicke theory with varying speed of light 
can solve the quintessence problem and it is possible to have a non-decelerated expansion of the 
present universe with BD-theory for anisotropic models without any matter. 
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^ I. INTRODUCTION 

Recently, the matter field can give rise to an accelerated expansion [1] for the universe 
stems from the observational data regarding the luminosity-redshift relation of type la 
supernovae [2,3] up to about z ~ 1. This matter field is called 'quintessence matter' (shortly 
Q-matter). This Q-matter can behave like a cosmological constant [4,5] by combining 
positive energy density and negative pressure. So there must be this Q-matter either 
■ neglected or unknown responsible for this accelerated universe. At the present epoch, a lot 

' of works has been done to solve this quintessence problem and most popular candidates 

for Q-matter has so far been a scalar field having a potential which generates a sufficient 
negative pressure. Since in a variety of inflationary models scalar fields have been used 
in describing the transition from the quasi-exponential expansion of the early universe 
to a power-law expansion, it is natural to try to understand the present acceleration of 
the universe, which has an exponential behaviour too, by constructing models where the 
matter responsible for such behaviour is also represented by a scalar field [6]. However, 
now we deal with the opposite task, i.e., we would like to describe the transition from 
I ' a universe filled with dust-like matter to an exponentially expansion universe and scalar 

fields arc not the only possibility but there are (of course) alternatives. Different forms for 
the quintessence energy have been proposed. They include a cosmological constant (or, 
more generally a variable cosmological term) a scalar field [4,7] a frustrated network of 
non-abelian cosmic strings and a frustrated network of domain wall [8,9]. All these proposal 
assume the Q-matter behaves as a perfect fluid with a linear barotropic equation of state 
and so some effort has been invested in determining its adiabatic index at the present epoch 
[10,11]. Recently, Chimento et al [12] showed that a combination of dissipative effects such 
as a bulk viscous stress and a quintessence scalar field gives an accelerated expansion for an 
open universe ( k = —1). Very recently, Banerjee et al [13] also have shown that it is possi- 
ble to have an accelerated universe with BD-theory in Friedmann model without any matter. 
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The possibility that the speed of light c might vary has recently attracted considerable 
attention [14-20]. In a cosmological setting, the variations in c have been shown to solve the 
cosmological puzzles - the horizon, flatness and Lambda problems of big-bang cosmology. 
The variations of velocity of light can also solve the quasi-lambda problems. For power-law 
variations in the velocity of light with the cosmological scale factors, Barrow et al [16] 
have shown that flatness problem can be solved. The Machian VSL scenario in which 
c = coa n , introduced by Barrow [15] has significant advantages to the phase transition 
scenario in which the speed of light changes suddenly from cT to c + , preferred by Albrecht 
and Magueijo [14]. For changing c, the geometry of the universe is not affected . We 
have allowed a changing c to do the job normally done by "superluminal expansion" . The 
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basic assumpssion is that a variable c does not induce corrections to curvature in the 
cosmological frame and that Einstein's equations, relating curvature to stress energy are 
still valid. The rationale behind this postulate is that c changes in the local Lorentzian 
frames associated with cosmological expansion. The effect is a special relativistic effect, not 
a gravitational effect. Therefore curvature should not feel a changing c. In a cosmological 
setting the postulate proposed implies that Brans-Dicke field equations remains valid even 
when c ^ 0. Magueijo et al [21] find the metrices and variations in c associated with the 
counterpart of black holes, the outside of star and steller collapse . The variation upon the 
theme are VSL theories which explicitcly break local Lorentz invariance, such as the one 
proposed by Albrecht and Magueijo [14] and for which black hole solutions remain exclusive. 

This paper investigates the possibility of obtaining a non-decelerating expansion (q < 0) 
for the universe in BD theory with varying speed of light in anisotropic models of the 
universe. 

II. FIELD EQUATIONS AND SOLUTIONS 

We consider the line-element of anisotropic space-time model 

ds 2 = -c 2 dt 2 + a 2 dx 2 + b 2 dfl 2 k (1) 
where a, b are functions of time only and 



dnl 



dy 2 + dz 2 , when k = (Bianchi I model) 

d6 2 + sin 2 9d(f) 2 , when k = +1 (Kantowaski-Sachs model) 

d9 2 + sinh 2 9d(f> 2 , when k = -1 (Bianchi III model) 



Here k is the spatial curvature index, so that the above three types [22] models are 
Euclidean, closed and semi-closed respectively. 

Now, the BD-field equations with varying speed of light are 



a b 8tv 
- + 2- = 

a b (3 + 2uj)<p 



(2 + w)p f +3(l+o;) % 



(2) 



db 8irp f 



k^_ 



a 
a 



uj I <fi 



(3) 



and the wave equation for the BD scalar field <p is 



8tt 



3 + 2cj 



Pf 



3p f 



(4) 



Here the velocity of light c is an arbitrary function of time, u> is the BD coupling parameter. 
p f and p f are density and hydrostatic pressure respectively of the fluid distribution with 
barotropic equation of state 

P f = (7 f - l)p f 

( j f being the constant adiabatic index of the fluid with < 7 / < 2 ). 
From the above field equations, we have the 'non-conservation' equation 
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As at the present epoch the universe is filled with cold matter (dust) with negligible 
pressure, so using p f = 0, the equation (4) and (5) gives 

and 

Now eliminating p f between (6) and (7), we have 

Since the scale factors, scalar field and velocity of light are time dependent, so we shall 
consider the following cases to obtain exact analytic form for the variables assuming some 
of them in polynomial form or in exponential form. 

Case I : The power law form of scale factors a, b and BD scalar field are assumed as 

a(t) = a t a ,b(t) = b/ 3 ^(t)=^ (9) 

where a 07 b ,(j) are positive constants and a,(3,p are real constants with the restriction 
a + 2(3 > 3 (for accelerating universe). 

If we use the relations (9) in (8), we have 

c = c t' 3 - 1 (10) 
where the positive constant c is restricted as 

2kc 2 (3 — 1) 

6g/i(a + 2/? + /i-l)(a + 2/? + /x-2)= (11) 
Also from (6) the expression for p f is 

P, = Po^- 2 (12) 

with 

Po = 3 gJ^V oMa + 2/3 + /i - 1) 
For consistency of the field equations the restriction between the parameters is 

(a + 2(3 + p £) 2 = M« + 2/3 + /x-l) [(2 + _ p) _ {a _ m)(3 + 2uj)] + 1 (ig) 



4 



Case II : The exponential form of scale factors and BD scalar field are assumed as 



a{t)=a e a \b{t) = b e l3 \<j>{t) = 4> e» t 
where a ,b ,<p are positive constants and a,(3,n are real constants. 

From (8), we have 

c = c„ ( 

where the positive constant c will satisfy 



(14) 



(15) 



(3 + 2uj)n(a + 2(3 + pf 



Ak(3c 2 



From (6), the expression for p f is 



p, = p e 



fit 



with Po = ^M« + 2/}^). 

From field equations the relation between the parameters becomes 



or (a + 2(3 + fi) <M + 



either a + 2(3 + fi = 
(3 + 2uj)/j,\ 



(P 



M(4 + 3w) 



(16) 



(17) 



(18) 
(19) 



In both the cases we have the deceleration parameter 

q = -l<0 

Case III : In this case , we have assumed the power law form of scale factor b, velocity of 
light c and BD scalar field <f> to be 



b(t) = b ^,c(t)=c t s ^(t) = <f> f i 
Using (20), we have from (8), 

rl 2 9fikr 2 



dt 2 



(3 + 2w)/x62 



(20) 



(21) 



To solve the differential equation we assume 8 — (3. It is to be noted that the above 
differential equation can be solved for all k. But for k = ±1 the solutions are not consistent 
to the field equations. So we consider only k — 0. In this case the explicit solution is 



u, + 2 2(^ + 1) ^ 

a(t) = a t 3-+ 4 +a 1 t 



2(w+l) 

b(t) = b t^+T- 

2(w+l) 



c(t) = c £ 3 "+ 4 



*3u,+4 



(22) 



with a , a x as integration constants. 
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In Figs. 1-5, we have shown the variations of q over t for various values of ui and the values of 
parameters ao = a± = 1. We have taken uj — —1.8, —1.5, —1, —0.87, —0.75 respectively in Figs. 1-5. 
In Fig. 6, we have shown the variations of q over t and uj in the range < t < lOand —2.5 < uj < 0. 



The deceleration parameter has the expression 



3[(3w+2)(a + Ql t) 2 +aJ(3w+4)t 2 ] 
+ [(3w+2)(a +a 1 t)+a 1 (3w+4)t] 2 



which is finite for all t and uj except t 



_o ]L (3a)+2) we have the non-decelerated 

universe for —2 < uj < —2/3. Also we have the non-decelerated universe asymptotically 
except for —1 < uj < —1/2. The variation of q over time has been shown graphically for 
different values of uj in figures 1-5. 
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III. CONFORMAL TRANSFORMATION : FLATNESS PROBLEM 

One important aspect of this model is that it can solve the flatness problem as well. To 
see this we make a conformal transformation [23] as 

9^u = (t>9nv (23) 

which enables us to identify the energy contributions from different components of matter 
very clearly. 

In this section, we have developed the BD theory in Jordan frame and to introduce the 
Einstein frame, we make the following transformations : 

dr] = ^a,d = ^a,b= ^b,ip = lii(p,p f = <j)~ 2 p f ,p i> = 4>~ 2 p p = ^ 2 Pj ,p^ = <p~ 2 p^ 

(24) 

So the field equations (2) - (4) transformed to 

a" J" A ( 3p,\ (3 + 2w) ,,, 

_ + 2j = -47T U + -f)- K —^r (25) 



b'\ 2 „a'b' kc 2 n (3 + 2uj) ,., 



and 



4) sib 

where ' = 

The scalar field ip (massless) behaves like a 'stiff' perfect fluid with equation of state 

p,=p, = Sg (28) 

If the velocity of light is constant, then in Einstein frame total stress-energy tensor is 
conserved but there is an exchange of energy between the scalar field and normal matter 
according to the following equation 



V (- 3p, 



c 2 



(29) 



On the other hand, if the velocity of light varies then we have separate 'non-conservation' 
equations 



/ a' b'\ ( p f \ V' ( 3Pf\ kcc' 



and 
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Thus combining the two energy densities, we have from the above equations, the equation 
for the conservation for the total energy is 

ft + 3-yHp = (32) 

Here £f=|^+2^is the Hubble parameter in the Einstein frame and 7 is the net 
barotropic index defined as 

7fi = 7A+7A (33) 

where 

n = n, + % = ^ (34) 

is the dimensionless density parameter. 

From equations (26) and (32), we have the evolution equation for the density parameter 

as 



Q' = n(n-l)[ 7 ff +2( 7 -l)ifJ (35) 
where H = t and H = \ . 

This equation in O shows that fi = 1 is a possible solution of it and for stability of this 
solution, we have 



7 < 3 (36) 

Since the adiabatic indices do not change due to conformal transformation so we take 
7j, = 1 (since p f = 0) and 7^ = 2. Hence from (33) and (34), we have 

7 = -J- — 37 

f $ 



Now due to upper limit of 7 , we must have the inequality 



n,<4|fij (38) 

according as 7 is restricted by (36). 

From the field equation (26), the curvature parameter Cl k = —kc 2 / b 2 vanishes for the 
solution Q = 1 . So for BD-scalar field it is possible to have a stable solution corresponding 
to O = 1 and hence the flatness problem can be solved . 



IV. CONCLUDING REMARKS 



We have performed an extensive analysis of solutions to Brans-Dicke theories with 
varying speed of light. For the power-law forms and exponential forms of the cosmological 
scale factors and scalar field in case I and // respectively, we have the velocity of light in 
the same form and we identified the cases where the quintessence problem can be solved 
for some restrictions on parameters. In the case I, we always have the accelerated universe. 



8 



For figures 1 and 2, q decreases to a fixed negative value asymptotically for u = —1.8 
and ui = —1.5 . For figure 3, q is linear with time and always negative for lj = — 1 . For 
figure 4, q decreases and then there is a singularity and after that q increases. In figure 5, 
q increases to a fixed positive value asymptotically. From figure 6, we conclude that the 
decelerated parameter is negative whenever — 2 < u < — 2/3 . In this range q decreases till 

t = — an d tncn increases to a fixed positive value asymptotically. 

Along with providing a non-decelerating solution, it can solve the flatness problem also. 
In fact, it has been shown that f2 = 1 could be a stable solution in this model. 
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